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ABSTimCT 


SOilCAP (Stiff Circuit Analysis Program) is a FORTRAN IV, Version 
2,3, computer program v^ritten for the CDC~6400~6600 computer series 
and SCOPE 3.0 operating svstem. It provides the circuit analyst a 
tool for automatically computing tiie transient responses and fre- 
quency responses of large linear time invariant nettjorhs, both stiff 
and non-stiff. The circuit description and user's program input 
language is engineer-oriented, making simple the task of using the 
program. 

Three volumes of documentation are available for the STICAP 
program; a theoary manual, a user’s manual, and a system's programmers 
manual. Volume I describes tlie engineering theories underlying 
STICAP and gives further references to the literature. Volume II, 
the user's manual , explains user interaction vjith the program and 
gives results of typical circuit design applications. Volume III 
depicts the program structure from a system's programmers viewpoint 
and contains fla,'/ charts and other software documentation. 
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CIPP^ER I 


r.ixrcDT?j!?iorT 

This document contaxns materxel describing tlie internal program 
structure of tiie STICAP program. 2\n overviex^ of the corposite system 
structure is gdven. The overlay structure is described, and the 
routines composing each overlay are included, for each routine gen- 
erated under contract IT?.Sl~9434-25 . A description of the COSXIAP and 
Gear routines may be found in the references [13 and [?] belov7, uhic.i 
contain as much documentation as is available for these programs. 

The program is written using CDC F0RTR2\H 4, Version 2,3. It is 
machine compatib3e ^rith the CDC 6400-6600 computer series, and runs 
under the SCOPE 3,0 operating system. It is segraentod in overlays of 
7CK or less using the SCOPE OVERLAY capability. All I/O is accom- 
plished using standard I/O files. 1/0 files are equivalenced so that 
Pile 5 is used for input and File 6 for output. Vo othe^^f'files are 
used by tins program. 


OUIGINAi; PAGE IS 
OP POOE QUAIOT 


1, Pottle, C., "Program Cornap - FORTRAtl Computer Routine,” System 
Theory Croup, School of Electrical Engineering, Cornell University, 
Ithaca, n, Y., 1968. 

2. Gear, C, 'T. , Algorithm 407 - DIPSUE for Solution of Ordinary 
Differential Eouations [D2] , Communications of the ACII, loarch, 
1971, Vol. 14, Ifo, 3, 105-190. 
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OISELAY L'XRUCTURS 

The first component of the program overlay structure {see Figure 
1) consists of a root segment (a main overlay) containing a super- 
visor routine STICAP, a comonly used routine USEFCIT, and blocks of 
labeled and unlabeled comiron core storage, also used by nany routines o 
The supervisor routine obtains control at initial eKecution and co- 
ordinates the calling into core of the successive overlays. This as 
accomplished by means of subsequent primary overlays, of vhich each 
successive one replaces the preceding, v/xth the number called depen- 
dent upon the mode of operation selected by tlie user. 

The first primary overlay consists of a block of routines which 
read the input circuit description cards, scaling card, and mode card. 
Tlie circuit description cards, whose format has been altered from 
that of COKIJFAP, provide the user an easy means of describing the cir- 
cuit. The scaling card, optionally present, determines the way the 
circuit equations are to be scaled for numerical computation of a 
solution. The mode card determines the program mode, or method by 
whidi t!ie solution is to be obtained. 

The second primary overlay consists of the original program 
COBHAP, minus the data input processing routines, \;hose function has 
been accon^lished by the routines mentioned above. The routines of 
this overlay process the circuit description data previously read, 
to obtain the mathematical equations governing the circuit. If the 
solution IS to be obtained by the COENAP mode, the original CORMAP 
routines continue processing tlie circuit equations. 

In the event that either the Gear or llatrix mode is selected, 
the third primary overlay is called into core. Data compiled by 
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Figure 1 - Overlay Structure 


^Appears xn a previous overlay also 
**Gear's ALGORITHfl 407 Routiiies. 



CORNAP IS compacted into a more economical form by the routine COI'IPAK. 
The control cards and data input cards for the mode selected are pro- 
cessed, and the final primary overlay is called. 

The final overlay routines are determined by the node selection. 
In one case the routines of Gear are used for numerical integration? 
in the other the matrix technique is used to obtain the solution. 
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SZSIIEU 0\?SRVIE?7 

An overview/ of the STICA? system indicacing the overall logical 
flov/ of the program appears xn the followxng pages, r?he fxgures 
indicate the routines T7ithin each overlay? the functions performed; 
and the calling sequences, or flo\; of control. 
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FUNCTIONS OF TUT!! IITDIVIdUMj ROUTINES 

In thxs section a brief c>escription of the function performed 
by each contractor supplied routine is given, and the routines are 
categorized by the overlay in vhidi tney operate., re remark tliat 
occasionally a routine may oe used in several overlays? hence it nay 
appear more than once» For a more complete description of a routine, 
the next section may be consulted.. 

The name and function of tlis various routines vill no\7 oe 
given . 

I . Root Segmenr Routines 

1. USEPCN - Supplies the user a means of specifying the 
independent sources of the circurh, 

2. STICAP - Supervisor routine ^’hich coordinates tlie calling 
into core of the primary overlays o 

II. First Primar^?^ Overlay 

1. REED - Processes title card and scaling card; coordinates 
routines \fhxch read CORJJhP circuit description. 

2. INPUT “ Reads circuit description cards in a free format. 

3. NODE ~ Reads Node select card and processes it. 

4. 2\P2\CK - Packs into one uord u-ie first character of each 
vjord from an array of uords. 

5. ATODSC - Interprets free format floating point numjDers. 

6. NMFJ - Scan routine for breaking out a name from a data 
field of free forr- data. 



*Some routines not contractor supplied are also indicated. Such 
routines are Langley routines or else Gear siipplementary routines » 



o 


7 . HODS - Scan routine for breaking out an integer number 
from a data field of free form dauao 

8. DECEQV - Changes alphanumeric BCD characters into internal 
integer format on the CDC ~ 6600 computer. 

9. OUTPUT - Reads and processes cards indicating the circuit 
outputs specified by the user, 

10. 2FSCAL - Reads scaling cards and stores scaling parameters, 

III. Second Primary Overlay 

CORNAP ~ The body of routines composing uhe original program 
CORNM*, minus routines i^hich read into core the circuit dsscrip- 
tion cards. This overlay also contains the logical decisions 
required for carrying out the various print options, such as 
state equations, transfer functions, etc, 

IV. Third Primary Overlay 

1. ?jPACK - Previously described? first primary overlay, 

2. CROLGR “ Coordinates the subroutines X7hich read the 
source order? initial conditions, run controls card 
group, and end card for t!ie Gear mode. 

3. ICONSV “ Reads and processes the initial conditions 
card groups for either the Gear or ilatrix modes, 

4. SRCORD - Reads and processes the source order card 
groups for tiie Gear modes. 

5. RTOTCTL - Reads and processes run control cards for 
Gear node. 

6. ATODEC - Previously described? first primary overlay. 
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7. DECEOV “ 'Previously c’eccribec’ “ first prii'ary overlay, 

0, CROLiDJ - Coordinates che oubrocrines ^jhich rea5 tlie source 
dsf ini cion, initial conditions, run controls card grouns , 
and tlie EI‘‘D card for the latrir rode, 

S, SRCiyT - Interprets the source definition card groun for 
the Ilatrir rode. 

10, RtJilCRL “ Reads and processes run conrrols card groim for 
the Hatri:: node, 

11, CQfTiYK - Rearranges data tables cutout bv CORUTi\P circuit 
translation rovcines into a rore econoiracal forr. 

12, CROLGi'I * Controllina routine for turd overlay, 

V, Gear Prir'aiy Overlay 

JL, IlfTGR - Coordinates Gear's routine for nurerical 

integration and determines ^fien the stop condition is 
reached , 

2. dIFSU'5 - Integrates the circuit equations in tlie Geer node, 

3, lj|iTPI;TT " Prints output data, 

4, DIFFUiT >- Svalua.tes tne rignt rier-Ciers of the state variable 
differential ecuations of the circuit, given a specified 
tine point, 

5, PEDERV - Evaluates the Jacooiaii Ilatri'i of t,ie svsten 
equation for the circuit. 

6. DECOIIP “ EU decor position routine used in tanden ”itn SOLVE 
to solve a linear syster', 

7. SOL’^TT “ Second stage of a LU decoi^nosition linear equation 
solver pair of routines, used by DI'>^SID, 
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VI, Matrxx Primary Overlay 

1. IIATRIX - Calls matrxx mode secondary overlays. 

2. SETSVE ~ A controllxng routxne for matrxx mode, \7hxch also 
obtaxns a sparse matrxx code. 

3. VAMINV - Inverts a Vandermonde matrix. 

4. EIGNV - Obtains system eigenvalues using QR transformation. 

5. HESEN - Reduces a matrix to upper Hessenburg form. 

6. QRTRAN - Performs a single double-step of the QR trans- 
formation, avoiding complex algebra. 

7. INTCL - Controls the closed form integration performed by 
the matrix mode. 

8. SVECOS - Processes integration of COSINE inputs. 

9. SVESIN - Processes integration of SINE inputs. 

10. SVESTP - Processes integration of STEP inputs. 

11. INIVMi - Integrates for impulse driving functions or other 

driving functions with a specified initial value. 

12. ^5PROD - Forms a matrix product. 

13. SAT - Calculates the matrix exponential exp (AT) and 
multiplies into a vector. 

14. SCEVAL - Evaluates the Du vector, u the vector of independent 
source inputs, y = Cx + Du the requested output vector, at 

a specified time. 

15. TRANS - Evaluates transient closed form solution at a 
specified time point. 

16. STEADY - Evaluates steady state closed form solution at a 
specified time point. 

17. STORE - Stores symbolic steady state closed form solution. 
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18. OUTPNT - Controls prxntxng of matrxK mode outputs. 

19. CLSOL “ Prxnts symbolxc closed form solutxons. 

20. TABPNT - Prxnts tabular form of matrix mode solutxons. 

21. AT2iK3 - Combxnes a sxne and cosxne functxon to obtaxn a 
sxngle sxne functxon. 



CILajPTJIJ V 

IHDIVXDUi^L ROUTIi'3 B_,SCRIPTIOiJS ' ^ 

Subroutines uSEPCiJ 

5?iiig routine allots cue user -cae capeubilicy of defining the 
fona of tne input vector for state-inpuc and stare ourput ecuariono 
of the neruork, for tne Goar node: 

dsr _ - 

— = Ak + Bu 
dt 

X7uere 

7 IS che state vector of the circuit; 

u IS tlie vector if independent sources, 

y IS tlie vector of outputs requested by the user. 

The user must either furnish PORTRAl'f equations defining the indepen- 
dent sources, or else muse supply a routine for computing their 
values, at a given time. The manner in \;hich this may be accomplished 
is described in the user's manual. 

Tne net effect of an entry to USEFCW must be tlie follouing; 

Given an initial tine T and a vector array U, compute values of the 
independent sources at time T and store as the components of the 
vector U. The order in vjnich tiiese components are stored must corres - 
pond to tlie order specified by the user '^/hen assembling the source 
oraer cards group (see user's manual). 

Inputs s a value of T and a vector array U 

Outputs . Tne values of tne independent sources at time T, stored 
in the array U, in the order specified by tne user’s source order 
cards . 


ORIGINAi; PAGE IS 
OF POOR QUALITY 
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Program 

Thxs program reads t!ie title card and calls IiJPUT, OUTPUT, and 
ZFSCAL to read the elements, output, and scaling card groups. The 
card groups are signaled by the folloizing header cards.* 

*ELE-XWTS 

*OUTPUTS 

*SCALIIJG 

where the asterisks appear in column one. The element and outputs 
card groups must always appear. The scaling card group need be 
present only if the circuit is to be scaled. The outputs card 
group must not be given before the elements card group. 

Tins program then calls I10D3 to interpret the mode card. The mode 
card must appear afuer tlie above card groups . 

The outputs from this routine are descrioed in the comment cards 
preceding the stLjroutines I&TPUT, OUTPUT, and ZFSCAL. 

SUBEOUTINE INPUT 

This subroutine reads tne elements card group. All cards read 
until the next card with an asterisk in column 1 are considered 
element cards. 

The cards are of the form 
NAI3 K1 N2 VM,UE 

for passive elements or independent sources and 
NAIIE la 112 VALUE*NAiS2 

for dependent sources. FAii3,ni,N2/ and VALUE must be separated by at 
least one blank. The asrerislc bett/een VALUE and iIAilU2 need not be 
present. IIAilE is the user ntune for the circuit element. Tnis name 
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laust be no inore tlian 4 cnaracters long, Tne firsr letter of t.ae 
element name tells the type of elo ^ent as follows. 

First Letter Element Type 

V Voltctge Source 

I Current Source 

R Pesisror 

L Inductor 

C Capacitor 

K Coefficient of Coupling 

II Mutual Inductance or Capacitance 

HI and IT2 for non~mutual elements are the nxamoer of the nodes between 
which the element is connected. The circuit nodes should be numbered 
compactly from zero, although the failure to do so is non~ fatal. 

HI IS given a positive reference wirh respect to H2. For mutual 
elements HI and H2 are the name of the t’-ro elements involved. For 
dependent sources NAIiE2 is tne name of tho controlling element i7ith 
a V or an I added to che beginning of the name depending on ^/hetlier 
the source is voltage or current controlled. 

A feature not normally employed by the average user is the 
follov7ing: 

A capacitor may be forced into the -proper tree- by placing 
the ■'/ord rree after the description of the element. 

HATOS HI U2 VALUE TIREE 

An inductor may be forced into tne cotree by similarly placing 
the V7ord cotree after the element description 
NAiSD HI H2 VALUE COTREE 



OUTPUT PAPJ^JSTERS 


* 16 - 


LIST 


K 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


12 


J 

1 

2 

3 

4-5 


4-5 

N F 

CON 6<1H 


A 3 t?ii'^ensional natnx ij.iose I^J^K eler«ent tslls 
infonratxon J of the Ith eleirent of type K. 

Type 

Dopendent Voltage Source 
Independent Voltage Source 
Zero Valued Resistor or Inductor 
Tree Forced Capacitor 
Capacitor 
Resistor 
Inductor 

Cotree Forced Inductor 
Zero Valued Capacitor 
Independent Current Source 
Dependent Current Source 
flutual Element 
Information 
User Marne for Element 
Elemenc Value 

Zero (later used to define outouts) 

For non-mutual elements the nodes oet^ieen 
i7hich tne element is placed 

For mutual elements the names of the fa7o 
elements involved 

A vector wiiose Ith element is the number of elements 
of type I (defined as in LIST) . 

An array naving two columns ^Jhose Ith row is the name 
of the element ^Jliicn controls tac Ith dependent 
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source. Column 1 is for der'cndent voltage sources 
and column 2 is for dependent current sources. 

ICON An array constructed like CONBCH, but \7hose eleraents 

tell the type of control of the dependent sources^. 1 
for voltage controlled and 2 for current controlled. 

MODEG A vector x^hosa Ith element tells tlie mamber of ele- 

ments and sources connected to the Ith node. 

HEX Tne number of nodes to xrhich no elements are connected. 

NODES The largest numbered node given by the user. 


Subroutine; MODE ( I! JODE , IP HINT ) 

Description ; This subroutine reads the mode card to determine x/hat 
program mode is selected and what print options are chosen. 

Program execution will be terrainated if an error occurs in 
describing the program mode. An error in describing a print option 
IS non-fatal, but the corresponding print request will be ignored. 
Input Parameters ; node Card - Tiis card will be in the follox\ring 
form 

*MODE KAfrlE, OPTION 1, OPTION 2, OPTION 3 
The card is in free form, except for an asterisk in colxn:\n one. Tug 
absence of the asterisk is non-fatal. Iiom^ NAME is the narxe of the 
solution program and is one of the follox;ing 
GEAR 
CORiAP 
IIATEi::, 

Tnese names ra}^ be alxbreviated or missnelled as long as the first 
letter of the first x?ord is correct. 
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The OiJtioaG follo'^jxiig che narie of -che solufxon ^>rogran are 
prxnt requests . The are separated hy cornrias and are of the 
follo\/xng forras 

SOLUTlOi'I EQUATIOKS 
STATE EQUATIONS 
TRANSFER PITNCTIONS 

These names may be abbrovxated or mxsspelled, as long as the 
first two letters are correct. The order xn whxch tlie prxnt op-cxons 
appear xs xnsignifxcant, nor do all or any need be present. 

Output Paraneters i IMODE - An xnteger tellxng what solutxon 
program to use 

IMODE = 1 COIttJAP 

IMODE = 2 GEAR 

IMODE == 3 MATRIX 

IPRINT “ A vector tellxng v/hau prxnt optxons '7ere requested. 

IPRINT(I) xs 1 xf optxon I x/as requested and 0 xf xt was 
not. The elements of IPRINT relate the follox7xngs 
IPRINT (1) State equations 
IPRINT (2) Transfer functions 

IPRINT (3) Solutxon equaLxons 

Addxtxonal Subroutxnes Reguxred 
APACK 

Examples of Usage 
1. Card read; 

*GEAR MODE 


Tne sujroutxne »^ould return 



ID 


IM0DJ3 = 2 
IPaiNT(l) = 0 
I?'HiTT(2) = 0 
IPRIJiTO) = n 

2 , Card read 

*MATRIX, SOLUTIOn EQUATIONS, TRAZJSPER FUI'TCTIONS 
The stCoroutins x^ould return 
IMODE = 3 
IPniUTU) = 0 
IPRIMT(2) = 1 
IPPJHT(3) = 1 

3 . Card read 

CORl^AP, TRM7SPER FUilCTIOHS 

Here the asterrs^c rn column 1 %7as om-cred. Tne sithroucine '70uld 
return 

IMODE = 1 
IPRIl'IT(l) « 0 
IPRINT{2) = 1 
IPRIHT(3) = 0 
4 „ Card read 

’’C,ST EQiI,TR 

Here the naities of tne soiutxon pronran and the prxnt requects ^jere 
aborevxeted. Tae Bixbrovitxne ’’oul<- return 
IIIODJ = 1 
IPRINT(l) = 1 
I?RIIJT(2) = 1 
IP?vIdT(3) = 0 
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Sxibroutine ; £PAC:C(KLET, ARRAY, M^ORD) 

Purpose s To pach the fxrst character of an array of \7ords into one 

word o 

Input Parameters ; APRAY - An array of x/ords containing alpha- 
numeric BCD characters. 

PXjET “ The nximber of x;ords in ARRAY 

Output Parairerers ; NI70RD - A xTOrd containing the first character 

of the first word in ARRAY as its first 
character, the first character of the 
second x/ord in ARRAY as its second 
character, and so on xmtii ARRAY is 
exhausted. The remainder of HJORD is 
filled v^ith blanks. 

Macnine Depend ancy s This su]:>routine assxuxies SO bit xrords x?ith 6 

bit alpha-numeric characters stored from the left. 

Example of Usage ; ARRAY (1) = 3HABC 

ARRAY (2) - IHD 

ARRAY (3) = 4HEPGH 

CALL APACK{3,ARRAY,l4TIORD) 

Upon return from the subroutine call, M40RD would contain SHADE. 

Subroutine5Required : 


PxSfIPT 

RSHFTA 


syst^ subroutines; not STICAP roxxtines. 


Subroutine ; ATODEC (WUFFER, LSRT, LSTP, DEC, HERR) 

Description ; This subroutine xjill interpret a nximber v^hich x/as read 
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t 

from a card and stored in an array in alpha-numeric BCD characters. 

The ntmber may be in integer, decimal or exponential form. 

The numloers and its exponent are presumed positive if no sign is 
given, numbers in exponential form need not contain a decimal point. 
All blanks idthin the input field are ignored. 

Tne subroutine starts at some given point in the input field and 
interprets until the first character not part of the nimber is 
found. This does not necessarily mean a non-nuneric character. 

The location of this character and the number found is outputted by 
the sxibroutine. If no number is found before a non-numeric charac- 
ter, such as a digit or a + or - sign, the nimber being interpreted 
IS assTimed to be a coefficient of 1.0. If no ntimber is found before 
the end of the input field, a 0.0 is given by the subroutine as the 
ntimber. 

Limitations : The magnitude of the number found must be beti/een 

10"*290 3^0+290 contain no greater than 15 digits. 

Input Parameters ; 

lItJFFER - A vector containing the card read in 80A1 format. 

LSRT “ A starting location in NDFPER x;here interpreting is 
to begin. 

Output Parameters : 

DEC - The number found. 

LSTP - The location in HDFFER of the next non-blank character 
follox7ing the number found. If the rest of the card 
was blan]:, LSTP = 81 
WERR - Error Code 


ITERR = 0 


Do error nresent 
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NERR = 1 Ilum-her is more than 15 digits long 

NERR « 2 Puinher is not between 10“^®^ and in 

magnitude. 

Machine Dependancy : There is no machine dependanev within tlie 

subroutine itself, ho’-jever a function \7hich it calls, DECEQV, is 
machine dependent. 

Examples of Usage : 

READ 10, (ilUPFER(I) ,1-1, 80) 

10 FORMAT (80 Al) 

(1) Card contains bbb2.147bbb. . . 

CALL AT0DEC(HUFFER,1,LSTP, DEC,FSRP) 

After the call 

NUPFER remains unchanged 

LSTP contains 81 

DEC contains 2.147 

KERR contains 0 

(2) Card contains bbA-2SIN(4T)bb. . . 

CALL ATODEC (NUPPER, 5 , LSTP , DEC , KERR) 

After the call 

ITUFFSR remains unchanged 

LSTP contains 6 

DEC contains 2.0 

KERR contains 0 

(3) Card contains bb21SE-3EXP (~5:)bb. , . 

CALL ATODEC (KUPPER, 1 , LSTP , DEC ,KEPJl) 
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After the call 

NUITPER remains unchanged 
LSTP contains 9 

DEC contains .216 

NERR contains 0 

(4) Card contains ~SIN3Xbbb 

CALL AT0DEC(NUPFER,1, LSTP, DEC, IIEPR) 

After the call 

NUFFER remains unchanged 

LSTP contains 2 

DEC contains -1.0 

NERR contains 0 

(5) Card contains bb3.1415926535897932bb. , . 

CALL ATODEC (NUFFER, 1 , LSTP , DEC , I^RR) 

After the call 

NUFFER remains unchanged 
LSTP undefined 

DEC undefined 

NEPJl contains 1 

{ 6 ) Card contains bb- 2 . 1 4E- 5bh 2 IS IN ( X) bbb . . - 

CALL ATODEC (NUFFER,!, LSTP, DEC, NERR) 

After the call 

NUFFER remains unchanged 

LSTP undefined 

DEC contains -2.14 


NERR 


contains 2 



24 


Subroutines reqtiirec! ; 

NCODE, DECEOV 

Subroutine : 

nai4E (MUFFER,LSRT,LSTR41SRR) 

Tnis stibroutine scans a card linage to find an element neime. 
This name must be less than 4 characters long. 

Input Parameters ; 

NUPPER The card image in 80A1 format. 

LSRT Starting location vrhere scanning is to begin 

Output Parameters ; 

I'lT'IORD Tne element name found. 

LSTP Location of blank after the element name. 

NSRR Error code . 

ttERR = 0 Ho error present. 

HERR - 1 !Iane over 4 characters long. 

HERR ~ 2 Ho name present (rest of card blank) . 

Additional Subroutines Reauired* 

— - --- - - - - — - — ---- - " -- 

APACK 
Subroutine ; 

NODE (NUPFER,LSRT,LSTP,HOM,NERR) 

This subroutine scans a card image from a given starting 
location to find the node number. The subroutine requires this 
node namloer to be followed by a blank. The largest node number 
must be less than 64. 
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Input Parameters ° 

NUFPER An array containxng the card xmage in 80A1 format. 
LSRT The starting location v/here scanning is to begin. 

Output Parameters; 

NUM The node number found. 

LSTP The location of the first blank follouing the 

node number. 

NERR Error code 

KERR ~ 0 Ko error present 

KERR = 1 A blank does not follo\/ the number. 

NERR = 2 Node number greater than 64. 

NERR = 3 Card after LSRT blank. 
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F\mction : DECEQV (NCOEP, LaTO) 

F;mction Purpose ; Convert an integer i7hose digits are BCD characters 
to an equivalent floating ooint number. 

Liiratations ; The integer being converted must be less than 15 
digits long . 

Input Parameters ; 

NCOEP - An integer array \7hose Ith element is the Ith digit 
of the integer being converted, read from the left 
in A format. 

LCNT ^ The number of digits in the number being converted. 
Output Parameters - 

DECEQV ~ The equivalent floating point number. 

Ilachine Dependency ; This function assmes 60 bit uords with alpha- 
numeric characters in CDC 6000 series BCD code stored on the left. 
Example of Dsage t 

HCOEF (1) « 1H2 
KCOBP (2) = 1H8 
IJCOEP (3) « 1H4 
NCOEP (4) = 1H5 
DEC = DECEQV (IlCOEP, 4) 

Upon return from the subroutine call DEC \^ould contain 2845.0. 
Subroutines Required ; 

ESHFT (a sysuem subroutine; not a STICAP routine) 

Subroutine ; OUTPUT 

This subroutine reads the output card group to determine the 
element currents and voltages selected as circuit outputs. This 
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siibroiitine t/ili te3aninate circuit processing if no outputs are 
selected or if all outputs contain errors. 

All cards read until a card \7ith an asterisk in colxnnn 1 is 
found are considered output cards. The cards are of the form 
VNAE1E 
or 

INA&3E 

depending on x^hether the voltage across or the cxirrent through WAI-JE 
is desired. Here NAME is the name of a circuit element given in the 
elements card group. 

Input Parameters ; 

HP - A vector of length 12 xvhose elements tell the number 
of each of the 12 types of elements contained in the 
circuit. 

LIST - A matrix x/hose (I,1,J) element tells the name of the 
Ith element of type J . 

Output Parameters ; 

LIST - A matrix x/hose (I,4,J) element gives the output informa- 
tion for the Ith element of type J, 

LIST (i,3,J)ssO No output taken from this element, 
LIST (I,3,J)al Voltage across this element is an 

output. 

LIST {I,3,J)='2 Current through this element is an 

output . 

Sxabroutine i ZPSCAL 

This subroutine reads the scaling card group. The circuit V7ill 
be scaled arovind a frequency or impedance value if one of the 
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follovrxng cards are present 
FREQUENCY = A 
liiPEDflNCE = A 

V7here A is a floating point number in free format. Scaling requests 
containing errors are ignored. 

Input Parameters ; None 
Output Parameters : 

Z SCALE - The impedance around t7hich the circuit is to be scaled. 
FSCALE - The frequency around T;hich the circuit is to be scaled. 
NUPPER - A vector containing the next header card. 

Subroutine : COtlPAK (EPS ,KA,KSC, NDERIV,NRONA, NLC/IIIODE, IP RINT, FSCALE, 

ZSCALE ,NAIILC) 

Description : This subroutine takes information assembled by CORNAP 

subroutines and puts it in the form used by the GEAR and IIATRIK mode 
subroutines . 

This siabroutine will terminate execution x;henever derivatives 
of the inputs are required by the output equations or v;henever deriva- 
tives of the inputs occur higher tlian the first degree, 

^Thenever the only derivatives of the inputs appear in the state 
variable equations and -then only the first degree, the state variable 
and output equations v;ill be of the form 

X = AX + Bu + 

Y = CX + Du 

The following substitutions are then made 

A • B1 + B -V B 


C • B1 + D -»■ D 
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The new circuxt equations no\j take the form 

^ = Aq^ + Bu 
y = Cg + Du. 

The names of the circuit outputs are augriiented by adding an I or V 
in front of the old name depending on \7hether the current through 
this element or the voltage across this element is the output. 

Inputs and Outputs ; All inputs are brought in by the t\JO arrays 
BLOCK and NAt'JES stored in common and by the argument list of the 
subroutine. 

The following variables are inputs to the subroutine but are 
not outputted, 

KA,KSC,NE0T7A *• Constants used to locate the h, B, Bl, C, D, 
and D1 matrices in BLOCK 

D1 - Coefficient matriK of source derivatives in the 

output equations 

The storage location of the remaining input and output variables are 
given in Figtrce 1. 

Machine Dependency ; The process of adding an I or V in front of the 
output names is machine dependent. 

The addition of the I or V is a means of more readily identifying 
current and voltage variables for printout and/or scaling purposes. 
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Subroutine s CR0L6R 

Description : This subroutine calls the stibroutines to read the 

source ordering, initial conditions, and nan controls card groups 
and tlie end card. These card groups are signaled by header cards 
which contain asterislcs in their first coluitin. The names on tliese 
header cards and the expected order in which they appear are: 

SOURCE ORDERING 

* INITIiil. CONDITIONS 

* RU1<I CONTROLS 

* END 

The source ordering, initial conditions, and run controls card groups 
need not appear in this order but must appear before the end card. 

In this order, hoirever, they t/ill process faster. 

The source order card group must appear if more than one source 
IS defined in the circuit description. Failure to do this X7ill re- 
sult in termination of program execution. 

If tlie run controls card group does not appear, the following 
run controls are assumed: 

ISTOP 1 
TO = 0.0 
TSTART =0.0 
H = 10“4 
HIIIN = 0.0 
METHOD = 1 
imiAX = 1 
MAXPTS =100 


INCRIN? = 3 
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A blank header card x;ill be assumed to be the end card. If the 
end card does not appear, the end of file card is assumed to be the 
end card. 

An unrecognizable header card t^ill result in termination of 
program execution. 

Usage ; Ihis subroutine is called by the main program x/hen the Gear 
option has been chosen, immediately after control is returned from 
COHNAP. 

Routine Uame ° ICOHSV 

Function ; Read cards containing the names and initial values of the 
coit^onents of the initial state vector. Print all nonzero conditions 
used by the program and/or any diagnostic error messages generated. 
Setting ; At execution time the ’* INITIAL CONDITIONS” card has just 
been processed, and the routine expects to read and process cards 
containing the initial circuit conditions, stopping at the first card 
read x/hich contains an * in column 1. 

Inputs ; 

Initial condition cards 

NSTV - The number of state variables COBKAP obtains 
NPSTV - The number of excess state variables* 

NA!1SV -- Names of CORHAP obtained state variables 
NAKLC - Names of excess state variables 
Outputs ; Error diagnostics and/or state vector of initial conditions. 
Components ii£ 5 )roperly foarmatted are set to zero. Excess components 
are ignored. 

*An excess state variable is a capacitor voltage or inductor current 
X'i'hose value is dependent upon other capacitor voltages or inductor 
currents . 
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Subroutine ; SRCORD 

Function t This subroutine reads the source order cards to determine 
a correspondence betireen the \7ay the sources are ordered in the state 
variable equations and the order in uhich they are computed and 
stored by DSEPQ'f. 

Inputs : 

WIN - The number of independent sources defined in the circuit 
description. 

NAIIES - A vector containing the user names of the independent 
sources in the same order they appear in the state 
variable equations. 

Outputs ; 

NORDER ~ A vector \7hose Ith element is the position in the 

Vector U(I) of USEFCP X7hich corresponds to the position 
of the Ith source in the state variable equation as 
ordered by the circuit translation routines. 

NUPPER •* A vector containing the next control card. 

Error messages are printed and the program terminated for the follox7- 
ing conditions: 

(1) An input is specified more than once. 

(2) An input specified is not found in "NMIES,” 

(3) Not all inputs defined in the circuit description are 
specified on the source order cards. 

Additional inputs : This s\fo routine X7ill read the source order cards. 

These cards contain the names of the sources, one per card in free 
format, in the order they are defined in USEFCN. All cards read are 
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assiimed to be source order cards, unless they contain an asterisk 
in column 1, t'Then a card containing an asterisk in column 1 is read. 
It IS stored in NUFFER and the source order cards are assumed to be 
conplete. 

Subroutine t RUWCTL 

Description . This s^lbroutlne reads the run controls card group for 
the Gear mode. The completion of this card group is signaled by a 
card uith an asterisk in column 1 or tlie end of file card. 

Input Parameters : The run controls cards nust be one of tlie 

following; 

INITIAL TII'IE = A 
STEP RESPONSE 
IMPULSE RESPONSE 
PRINT START = A 
STOP TIJ\E = A 
POINTS PRINTED = M 
EPS = A 

ADAT'^IS INTEGRATION 
HMAX = A 
HIIIN = A 
HINIT = A 

Here A is a floating point number and N is an integer. A description 
of the information conveyed by such cards is contained in the user’s 
manual for STICAP. 

Output Parameters ; 


TO - Initial value of time 
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TSTARP ~ Tiiae at which output printings is to begin 
ISTOP - stop condition 

ISTOP » 1 Plaximum number of points printed 
ISTOP = 2 Ilaximum value of time 
IIAXPTS - fiaximum number of points printed 
TSTOP - Ilaximum value of time 
H - Initial step size 
EPS - Error control 
I®THOD - Integration method 

IIETHOD = 0 Adams metliod 
METHOD = 1 Stiff integration metliod 
Ib^CRMT - Number of integration steps betv?een output printings 
IIAXDER “ Maximum derivative used by the integration method. 

IIAXDER IS set to 7 for lETIIOD « 1 and 8 for METHOD = 0. 
HIIAX - Maximum integration step size allox/ed 
HMIN - iiinimum integration step size allowed 
NPUT - Type of input 

NPUT =s 0 User defined 
NPUT = 1 Impulse 
NPUT = 2 Step 

NUFFER “ Vector containing the next header card. If an end of 
file card was encovintered , NUFFER (2) is set to IHE 
which will be interpreted as the end card. 

Default Options ; Unless specified by the rxm controls cards, the 
follo\7ing values are assumed. 

ISTOP = 1 
TO = 0.0 
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TSTART =s TO 
HPUT = 0 

HIIAX = (TSTOP-TO) /lO if TSTOP was specified, 

HMIN =0,0 

HilAX = 0,0 (interpreted as unbounded) 
imXPTS = 100 
EPS = 

METHOD = 1 
IWCRMT « 3 
H = 10~4 
imXDER = 7 

Machine Dependence ; This subroutine itself is machine independent 
but calls the machine dependent routines M>ACK and ATODEC, 
Subroutines Required ; APACK, ATODEC 

Program : CROLMX 

This subroutine calls the subroutines to read the source 
definition f initial conditions, and run controls card groups and 
the END card. These card groups are signaled by header cards which 
contain asterisks in their first columns. The names on these header 
cards are 

* SOURCE DEFINITIONS 

* INITIAL CONDITIONS 

* RUN CONTROLS 

* END 

The END card must be the last card appearing. If it is not present, 
the end of file card is interpreted as the END card. 
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The run controls card group nust be present. The source defi- 
nition card group must bo present unless the impulse or step response 
v/as specified in the run controls card group. If not execution t^ill 
be terminated. 

Subroutine ; SRCINT 
Description s 

This subroutine interprets the source definition card group for 
the matrix mode. This card group is signaled by a card x/itii an 
asterisk in column 1, follovjed by the words "SOURCE DEFINITIONS" in 
free format. The subroutine assumes all cards follov7ing this card 
up to the next card with an asterisk in column 1 are source defini- 
tion cards. These cards are of the form 
JJAI'IE = fl(t) + + ... + fN(fc) 

where NAME is the name of a source defined in tne circuit descrip- 
tion, and each fi(t) is one of the following 
A 

A*IMR 

A*SIN(B*T) 

A*SIK(B*T+C) 

A*COG(B*T) 

A*COS (B*T+C) 

where A, B, and C are floating point numbers in free format. The 
asterisks indicating multiplication need not appear. A source des- 
cription may be continued on additional cards by placing a dollar 
sign ($) in column 1 of the cards on \rhich the description is 
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continued as long as eacn fi(t) is completely described on only 1 
card and tfte total number of functions describing a given source is 
no more than 20. A sinusoidal or cosinusoidal function with a phase 
angle is broken into a sinusoidal plus a cosinusoidal function witn 
zero phase angle thus counting as two functions. 

This subroutine will terminate program execution if an error 
occurs in a source description or if not all the sources defined in 
tne circuit description are described in the source definition cards. 
Additional Input Parameters 

i-lIiSj - The number of circuit inputs. 

t^AMIN *• A vector containing the names of the circuit inputs 
in tno order they appear in the state variable equa- 
tions , 



NSORCE(K,J) = 3 cosine 


NSORCE(K,J) = 4 impulse 
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SOURCE - The magnxtude of cae Jtn function of the Ita source is 
contained in SOURCE (K,J, 1) / where K = NORDER(I) and 
the angular velocity {for sine and cosine functions) 

IS SOURCE (KxJ, 2) , 

Machine Oependance - This subroutine is machine independent but calls 
tne machine dependent subroutines APACK and ATODEC. 

Subroutine ; RtMCRL 

This subroutine reads in the run controls card group for uhe 
matrix method. 

The completion of this card group is signaled by a card i/ith an 
asterisk in colvimn 1 or the end of file card. 

Input Parameters 

The run controls cards read must be one of the following 
INITIAL Tirffi = A 
STEP RESPONSE 
IMPULSE RESPONSE 
PRINT START = A 
STOP TIME = A 
POINTS PRINTED = II 
PLOT I'NCP^iiENT « A 

v/here A is a floating point nximber and N is an integer. 

Outnut Parameters 



TO Initial value of time 

TSTART Time at which output printings is to begin. 

NPTS Number of output points to be printed. 

IJPUT Type of input 

DT Time between output printings. 
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Default Op cions 

If no xnxtial txme xs specxfxed, xt xs assumed to be zero. 

If no prxnt start time xs specified, it is assimed to be the 
initial time. 

If a stop time but no plot increment or points printed is 

specified or if only the plot increrient is specified, 100 
time points v/ill be printed. 

Machine Dependance 

This subroubine is inacnine independent but calls the machine 
dependent routines APACK and A.T0DDC, 

Program ; CROLGM 

This program calls COMPAK to rearrange the variables in blank 
common from the form supplied by CORNAP to the form required by 
GEAR and MATPJCX. One of the subroutines CROLGR or CR0LM:< is then 
called to read the card groups peculiar to the GEAR or the MATRIX 


mode. 
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Program IKTGR 

Thxs sabroutxne performs the xntegratxon of rhe state variable 
equations . 

Tae upper limit of the integration may be given by specifying 
the naximora value of the independent variable time or by specifying 
the number of points printed. 

If the maximum value of time is the stop condition, inter- 
polation IS used to determine the values of the state variables at 
the stop tine. 

Tne following error conditions will terminate program execu- 
tion 

1. Requested error is smaller than can be nandled for this 
problem. 

2. Corrector convergence could not be acnieved with H greater 
than HMIM. 

3. The requested error could not be achieved with K equal to 
Input Parameters 

IJSTV Tne number of state variable equations. 

TO The lo*v/er integration limit. 

XO A vector containing the values of tliG state 

variables measured at tirae TO. 

TSTART The time at v/aich printings is begun, 

TSTOP Tne upper Integra cion limit 

liAXPTS The maximum nujibcr of tine points to be printed. 

The stop condition 


ISTOP 
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HO 


ISTOP = 1 Maximum time points printed. 
ISTOP = 2 Maximum value of tine. 

The mrnber of circuit outputs 


HAMOIPT A vector containing tae names of uhe circuit outputs 

H Tne iniuial integration step size. 

HMIH The minimum integration step size. 

HMAX The maximiun integration step size permitted. 

IWCRMT The numJ^er of integration steps betveen printings . 

EPS The euclidean norm of a vector whose Ith element is 

uhe single step error of the Ith state variable divided 
by the maximum value of the Ith state varia.ble pre - 
viously encovintered must be less than this value. 

METHOD Tile integration metnod used. 

METHOD =0 An Adams predictor corrector method. 

METHOD = 1 A multistep method suitable for stiff systems. 

MAKDER The maximum derivative that should be used in the 

metnod. It must be less tiian 8 for Adams and 7 for 
stiff methods. 


Subroupine : DATPHT 

This si±>roUTL.ine prints the output for the Gear mode. It 
performs the calculation 

Y(t) = c x{t) + D u(t) 

^7here Y(t) is the output vector at time T. x(t) is a vector con- 
taining the values of the state variables at time T, and u(t) is a 
vector containing the values of the inputs at time T. 
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Input Parameters; 

T - present value of txne. 

X - a 8 X 30 matrix. x(l,I) contains the value of tne Itli 
state variable measured at tine T, 

C/D ~ constant matrices of che output equation. 

MO - number of circuit outputs. 

NSTV - number of stace variable equations. 

MIN - number of inputs defined in the circuit description. 

NORDER- a vector telling the order the inputs are defined in 
USEFCMo 

Output Parameters ° 

This subroutine prints the value of time and the values of the 
outputs measured at this time. 

Subroutine Required ; USEPCN 
llachme Dependency : None 

Suforouuine ; DI FFUM { T , X / DX) 

This Slab routine evaluates the derivatives of the dependent 
variables x, t7ith respect to tne independent variable T, using 
the matrix equation 
X=JxX+BU 

Input Parameters 

NSTV The numoer of stats variable equations. 

MIN The number of circuit inputs 

A Sy s tem matrix 

Input matrix 


B 
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T Present value of txne 

_X X(l/I) contains the value of the Itn state 

variable measured at time T. 

U U(I) IS the value of the NORDER(I)TH input 

measured at time T 

HORDER UORDER{I) is where the Kth source is located 

Output Parameters 

DY DY IS a vector containing X measured at the 

present values of X and T. 

Subroutines Required 
USEFa-I 

Soubroutine ; PEDERV 

Description ; This subroutine supplies DIFSUB i/ith the Jacobian 
ma-crix of the state variable equation; 

X = A X + B u 

which IS,, in fact/ the A matrix. 

Input Parameters ; 

M - nimoer of srate variable equations. 

A “ system matrix stored in core \/ith the first column of A 
in me 1st 30 storage locations, tne second column of A 
in the next 30 storage locations and so fortli. 

Output Parameters ; 

PW ~ Jacobian matrix. This is rhe system matrix. A, with 

each colurmi stored in core starting in the core location 
following the end of the previous column. 
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Subrou txnes Reauxred: Kone 

■ ■ ■i, m w, 

Machxne Dependancy ; None 
Example of Usage s 
Gxven tlie system matrxx 


A 


14 7 

2 5 3 : 

I 

3 6 9; 


A xs supplxed as 

27 blanks 27 blanks 810 blanks 

l/2/3/bfO/« • */b/4f5f6^b^b^ • • • /0^7 / Sf9/X)/bf o • * 

PW XS outputted as 


891 blanks 

1 , 2 , 3 / 4 , 5 / 6 , 7 , 8, 9 


Subroutxne s MATIKV{A,N, JSING) 

Purpose 

Invert a matrxx 
Descriptxon of parameters 

A - Input matrxx, destroyed xn computatxon and replaced by 
resultant xnverse. 

£I - Order of matrxx A 

JSING " Indxcates whether A xs sxngular; one x£ not, mxzius one 
xf xt xs- 

Remarks 

Matrxx A must be a general matrxx 



46 


Sxibroutines and fxmction subprograms reguxred: None 

Method 

The standard Gauss-Jordan method xs used. The determxnant 
xs also calculated. A determxnanfc of zero xndxcates that the matrix 
xs sxngular. 

Suaroutxne ; DECOMP (II, PS AVS^Jl) 

Function s 

This routine perforins an LU decomposition of the NxN matrix 
found in PSAVE. PSAVD is replaced by the matrices U and L, with U 
the main diagonal and above, L baloi’f the main diagonal. All main 
diagonal elements of L are unity, so need not be stored. 

If the PSAVE matrix is nonsingular, the routine returns a 
value of +1 for Jl? if singular, J1 is set to -1. 

The routines DECOMP and SOLVE are lifted directly from 

Forsythe and Moeller, Numerical Linear Algebra , p. 68 
and are well discussed in this reference. 

Subroutine ; SOLVE (K, PSAVE, SAVE, N5) 

Function . 

This routine employs the LU decomposition of DSCOtlP snored in 
PSAVE to solve for X the N equations 
(LU)x = b, 

with various b supplied by Gear's subroutine. The result is stored 
in SAVE. NS IS a parameter used for addressing the proper locations 
in SAVE in ^ihich x is to be stored 
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For large systems (K greater than 5) the LU decompositxon xs 
faster than the (^lATINV) matrxx xnversxon prevxously used by Gear 
If desxred, mxnor alteratxons to STICAP ^7x11 restore the Gear 
routine to xts previous form. These modifications are stated in 
Gear's ACM DIFSUB description (see STICAP theory manual for the 
complete reference) . 
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Program : tiATRlx 

Thxs overlay solves the system of state varxable equations by 
performing a spectral decomposition of the transition matrix then 
integrating in closed form. This technique is described in the papers 
"An Efficient llatrix Algorithm for the Solution of Linear 
Systems with X'^idely Separated Eigenmodes" ... by E. Young, 

J, Heii±>ocKel, and II. Ransom, Proceedings of the 4th Asilomar 
Conference on Circuit and System Theory, Pacific Grove, Calif., 
1970. 

Inputs to this program are described in the comments of the 
three secondary overlays and in the programmer’s manual. 

This first secondary overlay determines the eigenvalues of the 
system matrix, determines the inverse of the Vandermonde matrix given 
by these eigenvalues , and finally determines a sparse matrix code 
for the system matrix. 

CALL OVERLAY {3HCOI1, 4, 1,0) 

The next secondary overlay forms the transition matrix times 
a vector 

EXP (A*T) *V 

as a finite senes for the initial condition vector and for each 
driving vector. This expression is then integrated in close form 
and the symbolic solution evaluated at each time point at which out- 

put IS requested. 

CALL OVERLAY ( 3SCOH, 4,2,0) 

This final secondary overlay prints the tabulated results, and 
if requested, tlie close form solution. 

CALL OVERLAY (3HCOI1, 4, 3,0) 
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Program ; SETSVE 

*Ihxs program calls EIGW to deteinnine the eigenvalues of the 
system matrix and VAKINV to form the inverse of the Vandermonde matri 
given by these eigenvalues. Finally this program determines a sparse 
matrix code for the system matrix. 

Input Parameters ; 

A - The system matrix. 

N - The number of state equations. 

EPS “ The value belo^? which an element is felt to be zero. 

Output Parameters ; 

C - A complex matrix containing the inverse of the ’Jacobian 
matrix. 

KNZ - The number of non-zero elements in the system matrix. 
NZR,NZC - Vectors containing a sparse matrix code for the syster 
matrix. The Ith non-zero element of the system matrix 
is located in ro^>7 HZR(I) and colimn HZC(I). 

ANZ - A vector which replaces A in core v;hich contains only the 
non-zero elements of A. 

Subroutine ; VA 5 T 1 ETV 

This subroutine forms the inverse of the Vandermonde matrix 
given by the eigenvalues of the system matrix. The Vandermonde 
matrix of the n elements X3,/ \ 2 r ♦ • is the n x n matrix 


1 

1 

• • A X 

^1 


... Xn 



X 2 

• • • ''•n 



Q 

« 

• 

« 

• 

• 

• 

Xj^n-l 

X 

1 n-1 

... An 
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The xnverse Of the Vandermonde matrix is formed in this subroutine 
using a technique due to Kaufman.^ 


Let 


then 


CC3.3] = [Vi^l 

n-3 
I 

k=0 


-1 




*^13 * n 

n (Xi~Xk> 

k=l 

kj^i 

X 7 here the aj^ are the coefficients of the characteristic equation 

n 

P(X) s n (X-Xi) 

1=1 

In this subroutine the aj^ are formed from the recursive relationship 


Uq - L 

% = - i 2 : 


k Tn= 


^/here 


n 


m=l 


.ra 


^k-m *^m 


P=1 


He note finally that if “ ^ 31 ^ then = ^ 3 k** 

Input Parameters ; 

N - The mxnber of state variable equations. 

ROOT - A complex vector vjhose elements are the eigenvalues of 
the system matrix. 

Output Parameters ; 

C - A complex matrix containing the Vandermonde inverse. 


1. Kaufman, I., "Evaluation of an Analytical Function of a Companion 
Matrix with Distinct Eigenvalues," 1969 Proceedings of IEEE 
Letters, pp. 1180-1181, 
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Subroutine : HESEN(A,I-I) 

Supplier ; NASA Routine’^ 

Tins subroutine reduces a matrix through similarity transform- 
ations to an upper-Hessenfourg matrix. An upper-Hessenburg matrix u 
one in which all elements belox-? the subdiagonal are zero. 

This subroutine is the NASA Langley Research Center routine 
HESSEN. 

Input Parameters ; 

M - The order of the matrix being reduced. 

A - The matrix being reduced. This matrix is destroyed in the 
computations . 

Output Parameters ; 

A - The reduced matrix X7hich has replaced in input matrix in 
core. 

Subroutine ; QRTRAN(A,N,R,SIG,D) 

Supplier ; NASA Routine 

This stibroutine performs one <?ouble Q-R transformation on the 
matrix A. 

Let A = Ql^l 

^ ^-1 ~ Q 2®2 
^2 = "2Q2 = ”2'^ “ -2"^ 

\7here 

and Q 2 are unitary matrices 

and 

R^ and R 2 are right triangular matrices. 


*Not contractor supplied 
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A is then replaced xn core by A2. 

This subroutine is the NASA Langley Research Center routine 

QR2. 

Subroutine ; EIGNV 
Supplier : NASA Routine 

This subroutine finds the eigenvalues of a real, non-syiraaetric 
matrix using the double Q-R transformation. 

This routine was slightly modified from the NASA Langley Resear( 
Center routine REIS. 

Input Parameters : 

KJN - The number of state variable equations, 

AO - The matrix -whose eigenvalues are being determined. - 
Output Parameters : 

ROOT “ A complex vector containing the eigenvalues of the innut? 
matrix with real eigenvalues stored first and each com- 
plex eigenvalue followed immediately by its ccsnplex 
conjugate. 

Program ? INTCL 

This overlay calls the various integration routines to perform 

f 

the closed form integration for each function making up each circuit 
input. The close form solution is then evaluated at each time point 
requested. 

Input Parameters : 

W 


NPTS 


The number of state variables. 

The number of time points to be printed. 
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TSTART 

DT 

Af 8 /Cf D 

KIN 

NO 

XO 

INZ 

NZR,N2C,KNZ 

IPRINT 

NORDER 


ISORCE , SOURCE^ 
NSORCE 

NSC 

Output Parameters : 
OUT 


Kie time at which output printings are to be 

>< 

begun. 

The time increment between output printings, 
lllatrices in the input and output equations. 

The number of circuit inputs. 

The number of circuit outputs. 

A vector containing the initial conditions of ‘ 
the state variables. 

Tells \jhether all initial conditions are zero. 
INZ =» 1 if there are any non~zero initial 
conditions . 

Variables containing a sparse matrix code for 
the system matrix. 

Tells whether or not the closed form solution i 
to be printed. 

A vector relating the way the sources are 
defined in ISORCE, NSORCE and SOURCE V7ith the 
ijay they are defined in the output and state 
variable equations. 

Matrices containing information about sources. 
See SCEVAL for description. 

The number of different angular velocities 
exciting the circuit. 

A matrix containing the outputs evaluated 'at ^ 
each time point. 
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Subroutines required in thxs overlay ; 

STEJU5Y, INIVAL, STORE, TRMS, SCEVAL 
EAT, SVESTP, SVESIN, SVECOS , IIPROD 

Subroutine : SVECOS (AIIPL, 01 lEGA) 

TSiis subroutine determines the response of the system due to 
cosine driving functions. 

Input Parameters : 

N The number of state variable equations. 

NR The number of real eigenvalues of the system matrix. 

AMPL The amplitude of the driving function. 

OtIEGA The angular velocity of the cosine driving function. 

P A N3® complex matrix x;hose columns are the V’S of the 

expression 

EXP(A*T)*R s= V1*EXP(LAI©DA(1)*T) 

+ V2*F,XP(LAIIBDA(2)*T) 

+ ... 

+ VN*EXP(LAIIBDA(N)*T) 

where 

A IS the system matrix. 

R is the driving vector. 

LAMBDA (I) IS the Ith eigenvalue of A. 

and 

T 16 time. 

TO The time at v;hich the initial conditions of the system 

were measured. 

ROOT A conplex vector of length N containing the eigenvalues 
of the system matrix. 
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Output Parameters : 

CPCOS A NJQ? matrix containing the transient part of the 
solution . 

STYCOS A NX2 rectangular matrix containing the steady~state 
portion of the solution. 

Subroutine : SVESIN ( AI-IPL , OMEGA) 

This subroutine performs the integration of EXP (A{T-TAO) )*R 
*SIN (OMEGA*X) from TO to t uhere A is a matrix and R is a vector. 
Input Parameters : 

i 

i 

N The number of state variables. 

OMEGA The angular velocity of the driving function. 

R A vector of K length associated with the driving 

function. 

TO The initial value of time. 

ROOT A vector containing the circuit eigenvalues. 

P A matrix containing EXP(A*T)*R, 

NR The ntxnber of real eigenvalues. 

AMPL The amplitude of the driving function. 

Output Parameters : 

STYSIN A NX2 rectangular matrix containing the steady state 
solution to the integration. 

COEF A NXN matrix containing the transient part of the 
solution. 

Svbroutine : SVESTP (AMPL) 

This siabroutine determines the responses due to step driving 


sources . 
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Input Parameters ; 

N - The number of state variable equations. 

NR The number of real eigenvalues of the system matrix. 

ROOT A complex vector of length N containing the eigenvalues 
of the system matrix. 

TO The time at X7hich the initial values of the state variablei 
are measured. 

P A NXN complex matrix whose columns are the VS of the 
expression 

EXP(A*T)*R = V1*EXP(HD0T(1)*T) 

+ ... 

+ VN*EXP(ROOT(N)*T) 

where 

A is the system matrix 

R IS the driving vector 

ROOT (I) is the Ith eigenvalue of A 

and 

T IS time, 

Where complex conjugates occur, the complex conjugate V 
is omitted and the remaining VS stored consecutively in P 
AMPL The amplitude of the step function. 

Output Parameters ; 

COEP A NXN matrix containing coefficients corresponding to the 
transient part of the response. 

STP A vector of length N containing the steady state dc part 
of the response. 
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Subroutine : IHIVAL (AI5PL) 

This subroutine is called to determine the solution correspondir 
to the initial conditions or for impulse driving functions. 

Input Parameters ; 

N Number of state variables. 

ROOT Vector containing the eigenvalues of the system matrix, 

NR Number of real eigenvalues. 

P A matrix containing EXP (A(T-TO) ) where A is the system 
matrix. 

AIIPL The amplitude of the impulse function. 

Output Parameters : 

COEFO A matrix containing the closed form solution. 

Subroutine : I1PR0D(A,B,L,M,N) 

This subroutine multiplies a matrix A into a matrix B and 
stores the product in B. 

Input Parameters : 

A A matrix dimensioned 10x30. 

B A matrix dimensioned to have 30 rows but having a variable 
number of columns. This matrix is destroyed during the ' 
computation. 

L Number of rows in matrix A, 

M Number of colxmns in matrix A and rows in matrix B. 

N Nuiaber of columns in matrix B, 

Output Parameters ; 

B A matrix containing the product. 
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Subroutine : EAT (V) 

This subroutine finds the expression 


At 

e- 


V 


\jhere A is the system matrix and V xs a given vector. This expression 
13 given as 

At Xit X«t 

e~ • V = Vj^e ^ +. . .+ Vn® ” 


Here 


n 

Ya. « 2 

k=l 


« TT 


eigenvalues of A and C%j3 is the inverse 
of the Vandermonde matrix given by these eigenvalues. Since if 
« Xi*, = VjL* con^lex conjugate V’s are not given. 

Input Parameters ; 

K Tlie number of state variable equations. 

HR The ninnber of real eigenvalues in the system matrix. 

A A vector containing the non-zero elements of the system 

matrix. 

KHZ The number of non-zero elements in the system matrix, 

NZR,NZC T^ 7 o vectors containing a sparse matrix code for the 

system matrix. The Ith non-zero element in the system 
matrix is in row HZR(I) and column HZC(I) . 

C A complex matrix containing the inverse of the Vander- 

monde matrix given by the eigenvalues of the system 
matrix. 

Output Parameters ; 

p A matrix whose columns are the V's in tlie above expan- 

sion except that if Vj_ is complex/ Re[Vx3 is stored 
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in the Ith coluinn of P and Iin[Vi ] xs stored in the 

MU'*' 

I+lth column of P. ^i+1 which xs xs then omitted. 


Subroutine ; SCEVAL 

This subroutine evaluates the expression 


where D*U is the matrix expression in the output equation 
Y » C*X + D*U 

at each time point at which output is requested and adds the result 

to the contents of a matrix out which contains C*X at each output 

»• * 

time .point. 

Input Parameters : 

TSTART The time at which output printings are to start. 

DT The time increment beto^-een output printings, 

NUM The number of output time points to be printed, 

NIN The number of circuit inputs. 

NO The number of circuit outputs, 

NORDER A vector relating the way the sources are defined in 

ISORCE, NSORCE, and SOURCE VJith the way they are ordere 
in the state and output equations. 

ISORCE A vector whose Ith element tells the number of functior 
making up the NORDER (I) source. 

NSORCE A matrix vjhose 1,J element tells the type of the Jth 
function making up the NORDER (I) source. 

NSORCE ( I, J)=l STEP 

NSORCE ( I, J) *=2 SINE 

NSORCE (I, J)«3 COSINE 
NSORCE ( I, J)«4 inPULSE 
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SOURCE A 3 -dimensional matrix vhose element tells the 

magnitude of source HSORCE{I,J), and vjhose I/Jf2 element 
tells (for sine and cosine functions) the angular 
velocity of the function. 

D The matrix relating tlie outputs to the inputs in the 

output equation. 

Output Parameters ; 

OUT A matrix containing the circuit outputs at the requested 
time points. 

Subroutine ; TRANS 

This subroutine calculates the transient portion of the solution 
by stibstituting into the closed form expressions for the solutions. 
Input Parameters : 

N The number of state variable equations, 

NUM The number of different time values at which the output 

is to be measured, 

TSTART The time at which the first measurement is to be made. 

COEF A matrix containing the coefficients of the transient 
portion of the solution. 

ROOT A complex vector of N length containing the eigenvalues 
of the system matrix, 

NR The number of real eigenvalues of the system matrix. 

DT The interval between time samples. 

Output Parameters : 

OUT A matrix containing the outputs at the requested time 
points . 
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The values of the outputs are given at tirr»e 
TSTART+K*DT , K=0 , 1 , . . . ,tTuii-l 
Additional S\±iroutines Reouired; None 


Subroutine • STEM5Y (MTYPE , OMEGA) 

This sxibroutine calculates the steady state portion of the 

•V 

solution by stabstituting into tlie closed form solution. This cal- 
culation IS made at time 

TO + K*DT K = 0 , 1 ,... ,NUI1-1 

where TO is the time at which printing is to begin, 

MUM IS the number of time points at \7hich printing is 
to be made, 

and DT IS the printing increment. 

Input Parameters ; 

M The number of state equations. 

MUM The number of different time points at v/hicdi printing 
IS to be made. 

TSTART The time at which the first printing is to be made. 
NTYPE The type of driving function. 

NTYPE = 1 Sine or cosine driving function. 

NTYPE = 2 Step function. 

STYSC A matrix containing the closed form steady state portic 
of the solution for sine and cosrne driving functions.. 
STP A vector containing the closed form steady state portic 
of the solution for step functions. 


Output Parameters: 


OUT 


A matrix containing the circuit outputs at tlie request? 
time points. 
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Subroutine : STORE (OIIEGA) 

This subroutine stores the steady-state response for sine and 
cosine driving functions. The subroutine ^7ill coinbine responses witl 
the same angular velocities. 

Input Parameters ; 

HO The number of circuit outputs. 

STYSC A matrix \7hose 1,1 element contains the steady-state 

sine response for the Ith output and whose 1,2 element 
contains the corresponding cosine response. 

0MB(3A The angular velocity of the driving function. 

Output Parameters ; 

NSC The number of different angular velocities encountered 
so far. 

AHGVEL A vector i 7 hose Ith component contains the angular 

velocity of the Ith sine or cosine steady-state response 

SICO A matrix whose 1,J,K element contains the steady-state 
response of the Kth output ^;ith the Jth angular velocity 
and whose 2,J,K element contains the corresponding 
cosine response. 

Program : OUTPHT 

This overlay is used for printing output information for the 
matrix mode. Tabulated results are always printed. Close-form 
solutions are printed only if requested. 

Input Parameters ; 

IPRINT Tells whether or not the closed form solution is to be 
printed. 
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See CLSOL and TABPNT for additional input parameters. 

Subroutine ; CLSOL 

^ This subroutine prints the closed form solution for the 
requested outputs. 

Input Parameters ; 

NO The ntamber of circuit outputs. 

NA^OUT A vector containing the names of the circuit outputs. 

NSC Tne number of sine and cosine inputs. 

AHGVEL A vector x/hose Ith element is the angular velocity of 
the Ith sine or cosine input. 

SICO A matrix v/hose 1,J,I element tells the steady-state 
sinusoidal response of the Ith output due to the Jth 
sine or cosine input, and whose 2,J,I element tells tl 
corresponding steady state cosinusoidal response. 

STP A vector whose Ith element gives the steady-state stej 
response for the Ith output. 

N The number of state variables. 

NE The number of real eigenvalues, 

ROOT A complex vector containing the eigenvalues of the syj 
matrix \ 71 th the real eigenvalues stored first and eact 
complex eigenvalue folloxjed immediately by its comple: 
con3ugate. 

COEP A matrix containing the transient solution. The tran- 
sient solution for the Ith output is * 

COEP (I , J) *EXP tROOT ( J) *T] 

for ROOT{J) a real eigenvalue. fRien ROOT(J) is a comj 
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exgenvalue the solution corresoonding to ROOT(J) is com* 
bined vrith the solution corresponding to the complex 
conjugate of ROOT(J) to give: 

E}3>tRE(ROOT(J))*T] 

* [COEF(I,J) -SIN ( IM (ROOT (J))*T) 
+C0EP(I,J+1)*C0S(IM(R00T{J) )*T) 3 
v;here RE (ROOT (J)) is the real part of ROOT(J) and 
II1(R00T{J)) IS the imaginary part of ROOT(J) , 


Subroutine ; TABPNT 

This subroutine prints the tabulated results for the matrix 

mode. 

Input Parameters : 

TSTART Time at v/hich printing is to begin. 

DT Time increment bet\;een output printings. 

NRTS The number of time points to be printed. 

OUT A matrix whose I/J element is the Jth output measured 
at the Ith tire value. 

NO The number of circuib outputs. 

NAIiOUT A vector containing the names of the circuit outputs. 


Subroutine ; ATAIt3(X/Y) 

This subroutine comlDines a sine function X7ith a cosine function 
\7here both functions have equal angular velocities but no phase angl< 

a 

into a sine function with a phase angle. The phase angle is chosen 
to be positive or negative such as to mcake the magnitude of the phas< 
angle less than pi/2. 
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Input Parameters ; 

X The coeffxcient of the sxne functxon. 

Y The coeffxcxent of the cosxne functxon. 

> 

Output Parameters ; 

X The coeffxcxent of the resultxng sxne functxon - 
! 

Y The phase- angle xn radxans of the resultxng sxne functxon. 

Note: Both xnput parameters are destroyed xn the computatxons , 



CHAPTER VI 


FLOWCHARTS 

This chapter contains the flowcharts for all contractor supplied 
subroutines , 
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